Metastable states in the Blume-Emery- Griffiths spin glass model 
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We study the Blume-Emery-Griffiths spin glass model in presence of an attractive coupling 
between real replicas, and evaluate the effective potential as a function of the density overlap. We 
find that there is a region, above the first order transition of the model, where metastable states 
with a large density overlap exist. The line where these metastable states appear should correspond 
to a purely dynamical transition, with a breaking of ergodicity. Differently from what happens in 
p-spin glasses, in this model the dynamical transition would not be the precursor of a 1-step RSB 
transition, but (probably) of a full RSB transition. 



The Blume-Emery-Griffiths (BEG) spin glass model 1^,^ can be seen as a generalization of the Sherrington Kirk- 
patrick (SK) model, in which each site carries, beside to the spin variable Si — ±1, a lattice gas variable n-i = 0, 1. 
Its Hamiltonian is given by 

^ = - X! JijSiSjn.rij - ^ X! ^''^3 - l^^m (1) 

i<j i<j i 

where Jij are quenched Gaussian variables with zero mean and variance J'^ /N. The case K/ J = goes also under the 
name of Ghatak Sherrington model It has been found that in general, fixed the value K/J, there exists a tricritical 
point in the phase space. For higher temperature or chemical potential the model presents a second order transition 
from a paramagnetic to a spin glass phase. For lower temperature or chemical potential the transition becomes first 
order, with a jump in the mean density d — {ui) and overlap qap = {S^nfS^nf) between two replicas a and (3 (see 
Fig. li for the case K/ J = 0). Although the overlap parameter can become discontinuous, the transition seems to be 
quite different from the transition of other discontinuous spin glasses, like the p-spin models. A first obvious difference 
is that the transition of the BEG spin glass is first order also in the Ehrenfest sense, that is with a jump in the mean 
energy and entropy. Moreover, while in p-spin models the first step of replica symmetry breaking (RSB) turns out 
to be exact, in the BEG spin glass the question of the order of RSB appears much more subtle. The question was 
studied by many authors in the case of the Ghatak Sherrington model Q-H] . It was found that the replica symmetric 
solution is unstable against RSB in the whole spin glass region, with the appearance of complex eigenvalues in the 
stability matrix. Recently 0, it has been shown that this kind of instability persists also in the 1-step RSB solution. 
Therefore it seems reasonable that only a full replica symmetry breaking may probably give the correct solution, as 
it happens in the SK model. 

In the last years, a version of the model called frustrated Ising lattice gas (FILG), has been extensively studied as a 
model of a structural glass [p|~p^ . In this version of the model the disordered interactions are of the type Jij = ± J, and 
K = — J. Indeed, the FILG can be interpreted as a model for a system of asymmetrical molecules, with the variables 
Hi representing the presence of a molecule on the site i, and the variables Si representing the spatial orientation of 
the molecule. Two neighbor particles feel a repulsion of intensity 2 J if their relative orientation is "wrong" , that is 
JijSiSj < 0. Numerical simulations on finite dimensional lattices show that at low temperature the model develops a 
two step relaxation &jll|l, as observed in supercooled liquids and predicted by the mode coupling theory (MCT) of 
the glass transition ]16[|. 

Given these links between the model and the glass transition in structural liquids, it would be interesting to study 
(at least in mean field) if the discontinuous transition of the model is preceded at a higher temperature by a purely 
dynamical transition, from an ergodic to a non-ergodic phase, as it happens for p-spin glasses This would 

require the solution of the dynamical equations, as done for the p-spin, to see if they show a singularity, and of what 
kind, at some temperature greater than the static one. In this paper we address the problem from another point of 
view, exploiting the effective potential theory, introduced some time ago p9|-|2^ as a tool to identify the presence of 
metastable states, and a transition from an ergodic to a non-ergodic phase. We will refer to the so-called "annealed 
version" of the method. One considers two "real" replicas of the system (as opposed to the replicas generated by the 
replica trick), labeled 1 and 2, and defines a "degree of similarity", that is an overlap q, between them. Then one 
studies the system composed by the coupled replicas, with Hamiltonian H = Hi + H2 — e?- If F{e) is the free energy 
of the coupled system, then the Legendre transform V(q) = max[i^(e) + eq] will represent the "effective potential" 
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as a function of the overlap q, that is the work needed to keep the two repHcas, in the configuration space, at the 
distance q. What one finds in p-spin models [p^]-^ and in structural glasses in the HNC approximation p3| , ^ , 
and expects generally whenever metastable states are present in the system, that is for temperatures between the 
static and dynamical transition, is that the potential V{q) presents two minima. A lower minimum, for low values 
of the overlap q, corresponds to the "unbounded state" , in which the two replicas are in different "valleys" of the 
configuration space. A second higher minimum, for high values of q, corresponding to the "bounded state" , in which 
the two replicas are constrained to stay in the same valley. As the free energy of each single system is independent 
from which valley the system is in, the gap AV between the two minima can be interpreted as being equal to TS, 
where S is the logarithm of the number of valleys, and is called complexity or configurational entropy. 

We apply the effective potential theory to the model defined by ([l) with K/J = 0. The phase space is depicted in 
Fig. 1^, where the continuous curve represents the line of second order transition, the black dot the tricritical point, 
and the dashed curve the line of first order transitions. We now introduce two real replicas, coupled in the density 
overlap qd — jj X^i ^l^f- The Hamiltonian of the system is 

a— 1,2 \ i<j i j i 

where a = 1, 2 is the label of the real replica. By the usual replica trick, we can write the free energy of the system as 

= ^ (rf' + 9.') + ^ E + Z^) - ^ Tr exp(-Heff) (3) 

a</3 

where d — (n"") is the mean density, q^ = {S-^"- S'^°'n^°'n'^") is the spin overlap between two directly coupled replicas, 
QafS = {S°'"'S°'^n°'"n°'^) and Pap = {S^°'S'^^n^"n'^^) are respectively the overlap matrices between replicas with the 
same and with different "real label" , and TCcS is the effective single site Hamiltonian 

Heff = - (^d + (3^ E - f3'j' E ^-Z' E S^-S^f'n^-n^^ 

- 2/32 J2 E Pai3S^''S^^n'^n^^ - (3^J\, E " /^^ E W 

a<P a a 

In the replica symmetric approximation qa/s = q and pa/s = p. We make the further approximation q = p, and obtain 
for the free energy 

F = -2T\0g2+^{d'+q^,-2q^) 

-T Jvz log {1 + 2e" cosh(/3 J z) + e^^^^" [e^ cosh^ [PJs/q z) - sinh n] ) (5) 

where5= ^(d-g)+/3Af, n ^ (3"^ .P {q^ ~ q) , and 2?z g-^'/^. 

The physical states of the system will be given by the saddle points of the free energy, which are given by the 
equations 

r e" cosh(/3 z) + e^^+f^" cosh^ {PJJq z) - sinh Vl] , ^ 

d = Vz r ^ T — — T (6a) 

J 1 + 2e= cosh(/3J^z) + e2=+/'<^ [e*^ cosh^ (/JJ^gz) - sinhf7] 



e 



2H+/3£ cosh^ {(3. J yfqz) - cosh 17] 



T^z ^ ' / " 9 T (6b) 

1 + 2e= cosh(/3 J^g z) + e^'^+f^" [e^ cosh^ (/? Jy^ z) - sinh Q] 

^ f ^ / sinh(/3 jyg z) [1 + e5+/?^+» cosh(/3 J^g z)] \ ' 

* J ^ \^l + 2e"COsh(/3J7gz)+e2S+/3e[eJicosh^(/3jygz)-sinhn]y 

Notice that the saddle points are neither maxima nor minima of the free energy, because F has to be minimized with 
respect to d and qs, and maximized with respect to g. If e = then gs = g, 57 = 0, and one obtains the saddle point 
equations of the uncoupled system. Differentiating the free energy F with respect to e, we can evaluate the density 
overlap qd — {n^"n^"), 
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dF f e2=+'5'^ [e"cosh2(/3jygz) - sinhf]] 

~ ~ J ^ 1 + 2e^cosh(/3J^z) + e'^^+P^ [e^ cosh^ {(3 J ^ z) - sinhr^] 

In the above expression we can put e = = 0, and find the expression of the density overlap for a system without 
couphng, which cannot be obtained directly from the free energy of the uncoupled system. 

Now we want to compute the effective potential as a function of the density overlap qa- The procedure to do so was 
previously sketched, and consists in evaluating V{qd) — max[i^(e) + eqd]- We evaluate V{qd) in three points of the 

_ e 

phase space (marked by asterisks in Fig. [^) . One point at the edge of the region in which spin glass solutions of the 
saddle point equations exist, one on the line of first order transitions, and one in between. The result is shown in Fig. 

As soon as one enters in the region where spin glass solutions exist, a secondary minimum in the effective potential 
appears, signaling the existence of metastable states. We expect therefore that the dotted line of Fig. |l| will represent 
a line of purely dynamical transitions, below which ergodicity is broken and the relaxation function {ni{Q)ni{t)) will 
not decay to its equilibrium value, which in the paramagnetic phase is the density squared, but will stick at the value 
given by the secondary minimum of the potential. 

The presence of metastable states, and of a dynamical transition at a temperature higher than the static transition, 
in a model in which the low temperature spin glass phase is characterized by a full RSB, seems rather peculiar. It is 
reasonable to think that the usual picture, valid for p-spin models, according to which at the dynamical transition the 
phase space is splitted into an extensive number of valleys having all the same "distance" , that is the same mutual 
overlap, is modified in this case. Here in fact the spin glass phase is characterized by a complex ultrametric structure, 
typical of the SK model. In particular one may ask if in this case the gap /S.V in the effective potential can still be 
interpreted as a configurational entropy, that is the logarithm of the number of valleys in the configuration space. 
Finally, it would be interesting to understand how and to what extent this model can tell us something about real 
systems, that is structural glasses. 

I thank A. Coniglio, M. Nicodemi, M. Sellitto and J.J. Arenzon for many interesting discussions. This work was 
partially supported by the European TMR Network-Fractals (Contract No. FMRXCT980183), MURST-PRIN-2000 
and INFMPRA(HOP). 
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FIG. 1. Phase diagram of the model, for K/ J = 0. The region P is the paramagnetic phase, while the region SG is the spin 
glass phase. The full curve is the line of second order transitions, terminating at the tricritical point (black dot). The dashed 
curve is the line of first order transitions, and the dotted one is the line where non-paramagnetic solutions of the saddle point 
equations first appear. The three asterisks are the points where we have calculated the effective potential (see Fig. ^). 
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FIG. 2. Effective potential as a function of the density overlap qd, for chemical potential /i/J = —0.85 and temperatures 
T/J = 0.183, 0.119 and 0.0538 (points marked with an asterisk in Fig. |l|). The potentials for different temperatures are shifted 
so that the paramagnetic minimum is &t V = 0. The dotted lines join the potentials at the same temperature, and are only a 
guide to the eye. 
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